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1 Introduction 

All available experimental data on top quarks show that this particle behaves as a bare quark 
that does not hadronize. This fact implies that observables associated with or induced by 
top-quark spin effects can be used as reliable tools for investigating the interactions of this 
quark. As far as hadronic top-antitop quark (tt) production is concerned, the correlation of 
the t and t spins and polarization of the t and t samples in tt events were first experimentally 
explored at the Tevatron [1-4]. Evidence for ti spin correlations was first reported by the 
D0 experiment [3] and they were first observed by the ATLAS experiment [5]. The ATLAS 
and CMS experiments at the LHC (7 and 8 TeV) have performed a number of top spin 
correlation and polarization measurements, using various sets of observables [5-11]. These 
measurements have reached a level of precision that allow for a meaningful interpretation. 
All measurements of ti spin correlations and top-quark polarization so far are in agreement 
with respective Standard Model (SM) predictions [12-15]. 

The degree of top-quark polarization in a ti sample and the strength of the correlation 
of the t and i spins depends on the production dynamics and, for a ’given’ production 
dynamics, on the choice of observable, that amounts to a choice of reference axes. (These 
axes may be interpreted as t and i spin quantization axes.) In order that future measure¬ 
ments at the LHC will be able to extract complete information about the ti production 


- 1 - 




dynamics from their data as far as top-spin effects are concerned, it is desirable to devise 
a set of observables, in some basis, which probe all of the independent coefficient functions 
of the top-spin dependent parts of the tt production spin density matrices. In this paper 
we present such a set. The members of this set are constructed such that they have def¬ 
inite properties with respect to discrete symmetries. For dehniteness we consider, as far 
as tt spin correlation effects are concerned, the dileptonic tt decay modes and, as far as 
polarization effects are concerned, also the lepton plus jets hnal states. These modes are 
the most sensitive channels for the investigation of top-spin effects. We compute the distri¬ 
butions and expectation values of these observables at next-to-leading order (NLO) QCD 
with respect to ti production and decay and take into account also the NLO mixed QCD 
weak-interaction contributions and for the transverse t and t polarizations also the mixed 
QCD-QED corrections of order a'^a. In addition, we assess the sensitivity of our set of 
spin correlation and polarization observables to new physics (NP) effects in the framework 
of the effective held theory approach. We take into account all gauge invariant operators 
of dimension six which are relevant for hadronic tt production [16-22]. We compute the 
contributions of the associated anomalous couplings to our observables. There is already an 
extensive literature on phenomenological studies, including [23-48] and [15, 49-68, 70, 71], 
on how anomalous couplings, in particular, a non-zero chromo-magnetic and/or chromo¬ 
electric dipole moment of the top quark would affect ti production at hadron colliders. In 
fact, if one considers only the two chromo-moments and puts all other anomalous couplings 
to zero, it suffices to take only the experimental results for the ti cross section at the Teva- 
tron and the LHC (7 and 8 TeV) and the respective NNLO QCD predictions [72] in order to 
obtain sensitive bounds on these two anomalous couplings [50, 55, 68]. The point we want 
to make with this part of our analysis is that our set of observables, including the ti charge 
asymmetry^ is i) large enough to get information on all anomalous couplings of our effective 
Lagrangian (in fact the number of observables in this set exceeds the number of couplings) 
and ii) these observables allow to disentangle the various anomalous contributions to a 
considerable extent. 

This paper is organized as follows. In the next section we decompose the ti production 
density matrices of the reactions gg, qq —)• ti in terms of an orthonormal basis and classify 
the coefficient functions of these matrices with respect to parity P, CP, naive time reversal 
Tn, and Bose symmetry. The decomposition applies to all parton reactions ij —)• ttX. In 
section 3 the effective Lagrangian £np is introduced that we use for describing non-resonant 
new physics in hadronic ti production. We use the operator basis employed in [22], albeit 
with different normalization conventions for the anomalous couplings. Possible new physics 
contributions to the t ^ Wb decay vertex do not affect our leptonic angular correlations 
and distributions dehned in section 4. We briefly recall the presently available bounds on 
the anomalous couplings contained in Tnp and discuss the domain of applicability of the 
linear approximation used in section 4. In this section, we consider ti production and decay 
into dileptonic and lepton plus jets hnal states at the LHC. For the dileptonic hnal states, 
we present a set of leptonic angular correlation variables that are constructed such that each 

^The tt cross section should, of course, also be used but is, in principle, redundant; see section 4. 
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member of this set probes a different tt spin correlation coefficient of the production density 
matrices. In addition, we list a set of single lepton distributions/observables which probe 
the t and t polarization coefficients of the production density matrices. These observables 
can be applied both to dileptonic and lepton plus jets final states. We compute, for tt 
production at the LHC at 8, 13, and 14 TeV, the two- and one-dimensional distributions 
associated with these observables at NLO QCD (tt production and decay) including the 
NLO mixed QCD weak-interaction corrections, for the transverse t and i polarizations 
also the mixed QCD-QED corrections, and, in addition, the contributions induced by Tnp 
which are linear in the anomalous couplings. For completeness we determine also the 
contributions of £np to the ti cross section and the ti charge asymmetry at the LHC, 
Ac- For the observables considered in this paper we estimate the application range of 
the linear approximation in the anomalous couplings. Our results show that our set of 
correlation and polarisation observables, including the tt charge asymmetry Ac, allows 
to determine all anomalous couplings of Tnp and disentangles these contributions to a 
considerable extent. We conclude in section 5. 


2 tt production density matrices in an orthonormal basis 

We consider in this section, for the purpose of setting up our notation, hadronic tt produc¬ 
tion by the 2 —)> 2 reactions 

9{Pi) + a{P 2 ) -t t{ki,si) -k i{k 2 , S 2 ) , (2.1) 

aiPi) + q{P 2 ) -t t{ki, si) -k i{k 2 , S 2 ), (2.2) 

where Pj,kj and si,S2 refer to the 4-momenta of the partons and to the spin 4-vectors of 
the t and i quarks, respectively. The production density matrices of these reactions are 
defined by 

-^ 0102 , /3i/32 ~ E {tiki, a2),t{k2,l32)\T\aipi), 6 (^ 2 ))* 

X {t{ki, ai),i{k2, Pi)\T\a{pi),bip2)) (2.3) 


where I = ab = gg, qq and the bar denotes averaging over the spins and colors of I and sum¬ 
ming over the colors of t, i. Moreover, a, f3 are spin labels referring to t and i, respectively. 
The matrices can be decomposed in the spin spaces of t and t as follows: 


R' = fi 
fgg ~ 


H^ll (g) 11-k 4^+0-* 

(dTras)^ 

Nc{N^ - 1) ’ 


]l + 4^-]l0a* + CLa*0cjJ 


fgg — 


iV2 


(2.4) 


where W denotes the number of colors. The first (second) factor in the tensor products of 
the 2x2 unit matrix 11 and of the Pauli matrices cj* refers to the t (t) spin space. 

The functions and C-j can be further decomposed, using an orthonormal basis 
which we choose as follows. Here and in the following k denotes the top-quark direction 
of flight in the tt zero-momentum frame (ZMF). Here p = pi denotes the direction of 
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one of the incoming partons in this frame. Then the following set forms a right-handed 
orthonormal basis: 


{f, k, n} : f=-(p-?/k), n=-(pxk), y = k p, r = \/l - yT ( 2 . 5 ) 

Using rotational invariance we decompose the 3-vectors and the 3x3 matrices C-j 
(which have a symmetric and antisymmetric part with 6 and 3 entries, respectively) with 
respect^ to the basis (2.5): 

bI^ = bi^ h + h + bi^ hi , ( 2 . 6 ) 

(5/j = 4^ hfj + 4^ kkj + 4n ^ihj 

+ 4k 4ikj + kirj) + 4n {kihj + hikj) + 4n + nifj) 

+ ^iji{4ki 4k + 4'^i) ■ ( 2 - 7 ) 

The coefficients b^^, 4^/ are functions of the partonic c.m. energy squared, s, and of y = k-p 
which is equal to the top-quark scattering angle cos0j. The terms in the antisymmetric 
part of (2.7) may also be represented as follows: 

4^ijiri = 4{kinj - hikj), 4^ijiki = 4{nifj - nhj), 4^ijini = 4{nkj - hf-j). ( 2 . 8 ) 
Because the initial gg state is Bose symmetric, the matrix must satisfy 


R33{-p,k) = R39{p,k). 


(2.9) 


If CP invariance holds then 

^iia2,/3i,/32(P’'^) = ^^i,/32,«i«2(P’>^)’ i = sg,qq- ( 2 . 10 ) 

The implications of (2.9), (2.10) and of parity invariance on the structure functions defined 
in (2.6), (2.7) are given in Table 1. In this table the properties of these functions with respect 
to a naive Tn transformation (reversal of momenta and spins) and CPTn transformation are 
collected, too, where just for this purpose, the absorptive part of the respective scattering 
amplitude has been neglected. The entries in Table 1 are to be read as follows: For example, 
for a T-invariant interaction like QCD one has bl4{y) = —hl4{y) and thus b^k = 0 at Born 
level, whereas absorptive parts generated at 1-loop render this function non-zero. The above 
decomposition of the spin density matrices R^ and classification with respect to discrete 
symmetries is analogous to the analysis of [73] where a non-orthogonal basis was used. 

The spin density matrices of the 2—7-3 reactions that contribute to ti production at 
NLO in the gauge couplings and are inclusive in ti have the same structure as those discussed 
above (with coefficient functions which depend on more variables) and can also be classified 
with respect to discrete symmetries. Table 1 is useful as a guide for the construction of 
observables at the level of dileptonic and lepton plus jets final states that project onto 
contributions to the differential cross section which are even and odd with respect to the 

^The decomposition of the Clj is not unique because 5ij = hifij + fifj + kikj. 
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Table 1: Transformation properties of the structure functions defined in (2.6), (2.7) with 
respect to CP, P, and Tn- The implications of Bose symmetry in the last column apply 
only to / = gg. 



above discrete symmetries. These observables apply also to investigations beyond LO, cf. 
section 4. 

Strictly speaking, a classification of observables with respect to the discrete symmetries C 
and CP is not possible in proton proton collisions because the initial pp state is neither an 
eigenstate of C nor of CP. Yet, for the main tt production modes, gg, qq ^ tt + X it can be 
done. The ‘contamination’ of CP-odd observables by contributions from parity-odd mixed 
QCD weak contributions to gq{q) —)• tiq{q) is at the level of a few per mill [14]. 

We add a few remarks on the coefficient functions and on the information contained in 
Table 1. 

i) The functions and are generated by P- and CP-conserving 

interactions, in particular by QCD already at tree level. The functions determine 
the tt cross section, the pip distribution, etc., while the functions {a,b = r,k,n) 
induce P- and CP-even tt spin correlations. 

ii) The functions 6^, are generated by P-violating interactions and induce longitudinal 
t and i polarizations (i.e., a polarization in the production plane). Non-zero differences 
b^ — b~, b'^ — b'^, which are P-, CP-, and CPTx-odd, require CP violation and 
absorptive parts. 

iii) The coefficient functions b^ are generated by P- and CP-even absorptive parts in the 
scattering amplitudes, which is the case at 1-loop in the SM, and already at tree level 
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by P-invariant, but CP-violating (effective) interactions. If interactions of the latter 
type are present then — b~ 7 ^ 0. The functions b^ induce a transverse t and t 
polarization (i.e., a polarization normal to the production plane), which would differ 
for t and t if non-standard CP-violating interactions are present. 

iv) Non-zero coefficient functions Ckn, Cm require P-odd but CP-even absorptive contri¬ 
butions to the scattering amplitudes. In the SM they are generated by the absorptive 
parts of the 1-loop mixed QCD-weak corrections, and they are very small, see sec¬ 
tion 4. 

v) The coefficient functions Cr, Ck are generated by P- and CP-violating interactions. 
As is well known, in the SM CP violating effects in flavor-diagonal reactions like 
ab — 7 - tt generated by the Kobayashi-Maskawa phase are tiny. Observable CP effects 
in hadronic tt production requires sizeable NP CP violation. 

vi) A non-zero coefficient function Cn requires P-even but CP-odd absorptive contribu¬ 
tions to the scattering amplitudes. In the SM this coefficient function is tiny. 

As an aside we recall that in a time-reversal invariant theory, expectation values of Tn- 
odd observables can be generated by the absorptive part iAfi = Tfi — 71)- of the T-matrix 
element Tfi of a reaction i ^ f. T-invariance and the unitarity of the S matrix imply 

- \Tfi\^ = -2im(M/.r/j -, 

where i, / denote the initial and final state with 3-momenta and spins reversed. The left- 
hand side of this equation has the generic form of a T^-odd observable. 

3 Effective NP Lagrangian for hadronic top quark pair production 

We assume that new physics (NP) effects in hadronic tt production and decay are induced 
by new heavy particle exchanges, characterized by a mass scale A. Then one may construct 
a local effective Lagrangian Cnp in terms of the SM degrees of freedom which respects 
the SM (gauge) symmetries and describes possible non-resonant new physics interaction 
structures as long as the moduli of the kinematic invariants of the tt production and decay 
processes < A. Recent analyses in the context of tt production include [18-22]. Here, we 
follow [22], but take also CP-violating effective interactions into account. 

3.1 tt production 

We concentrate on new interactions described by operators of mass-dimension dimO < 6 
which manifest themselves in hadronic ti production. We consider only that set of operators 
which contribute to tt production through interference with the dominant SM amplitudes, 
i.e., the QCD tree-level amplitudes. 

Let’s first list the operators which involve gluon fields. We do not consider operators that 
affect only the pure gluon vertices, such as GGG or GGG, where G (G) denotes the (dual) 
gluon field strength tensor. For a recent bound on the anomalous coupling associated with 
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GGG, see [70]. The strength of GGG is severely constrained by the upper bound on the 
electric dipole moment of the neutron. For dimO < 6 one then remains with the following 
three gauge-invariant operators which involve t,t and up to three gluon fields [18, 20, 22]: 


Ogt = [tR^T^D’^tn] G%, (3.1) 

OgQ = (3.2) 

OCDM = G%, . (3.3) 

Here Ql = is the left-handed third generation doublet and $ = ^ 172 $^ = {4>o, —0-) 

is the charge-conjugate Higgs doublet field. Moreover, = dg, + iggT'^G^^ and = 
dg_G^ — dyG^^ — gsf^'^^G’^^Gy is the gluon field strength tensor. Furthermore, r“ are the 
generators of SU(3)c in the fundamental representation, with ii{T°‘T^) = bahl‘^- 
Using the equation of motion for the gluons one can show [22] that the sums Ogt -|- 
and OgQ -|- OgQ are given by linear combinations of 4-quark operators given below, and are 
thus redundant if one uses the 4-quark operators below as independent set. On the other 
hand, the combinations Ogt — Og^ and OgQ — G^gQ induce CP-violating interactions, which 
cannot be expressed in terms of the the 4-quark operators below. One can therefore use a 
hermitean effective Lagrangian with these two operator combinations and with (3.3) that 
contains three independent complex coefficients. After spontaneous symmetry breaking, 
(4>) = w/\/2, {v — 246 GeV) and restriction to the top-quark gluon sector, one obtains, using 
operators with definite P and CP properties, real and dimensionless coupling parameters, 
and the (on-shell) top-quark mass mt for setting the mass scale^: 


-^NP,g 

where 


and 


9s 

2mt 






m-t L 




Of__) = 0.+0l, Of_^^ =i{o+- oi), 


(3.4) 

(3.5) 


0+ = O- = G 


flU 


(3.6) 


Here t denotes the top-quark Dirac field. In (3.4) the real, dimensionless parameters jit and 
dt denote the anomalous chromo-magnetic (CMDM) and chromo-electric (CEDM) dipole 
moments of the top quark in the convention of [15]. Contrary to our previous analysis in 
[15] we stay here strictly within the effective field theory approach - i.e., we do not interpret 
these moments as form factors which may have an absorptive part^. The chromo-magnetic 

dipole operator is P- and CP-even, O^ _is P-even and CP-odd, while the chromo-electric 

dipole operator and O®_^ are both P- and CP-odd. 

Next, we list the set of gauge-invariant dimO = 6 four-quark operators that yield 
non-zero tree-level interference terms with the QCD amplitude for qq —)• tt. We assume 
universality of the new interactions with respect to the light quarks q ^ t. Because the 

®Our sign convention for the QCD quark gluon interaction is ^.9 = -gsq-i^T'^qGl. 

^In [22] only the chromo-magnetic dipole operator was considered. Their parameter Chg is related to jit 
by uRe Chg/{V2A'^) = -gsjit/{2mt). 
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respective parton distribution functions suppress tt production by qq annihilation for q = 
s, c, b as compared to q = u, d, we take into account only the contributions with u, d in the 
initial state. 

One is then left with the following 7 gauge invariant operators [20, 22]: 


[QLl^.T-QL){qLl^T-qL), 

(3.7) 

o^Qf= (QL7/.rv^gL)(gL7^rv^gi), 

(3.8) 

{iRl^.T<^tR){uRrT^UR), 

(3.9) 

Ou= {iRlt.T^tR){dRrT'^dR), 

(3.10) 

[QLlt.T-QL)[uR^^T‘^UR), 

(3.11) 

0g= {QLlf.T^QL){dR7>^T<^dR), 

(3.12) 

Ojf = {iRj^THR){qL7>^T\L). 

(3.13) 

{ul,ul) denotes the first generation quark doublet, and 

denote the Pauli 


matrices {tv{a^a^) = 25 ab-) 

These operators contribute to the effective NP Lagrangian a term of the form 
f8 ll 

denoting Oq^ = Oi, etc. Again, the terms which involve b quarks are not of interest to us 
here. As pointed out in [22] it is advantageous to combine these seven operators such that 
one gets four isospin-zero operators with definite P and C properties and three isospin-one 
operators. (In the isospin-one case it is not possible to combine the operators such that 
they have definite properties with repect to C and P.) In the following q = {u,d) denotes 
the isospin doublet. One gets 

-^^NP,q = -Cnp,o + -^^NP,1 , (3-14) 

where the isosoin-zero part is 

-^^NP,0 

and 

Oyy = (g7^r“g) , (3.16) 

Ova = (qi^T^q) , Oap = {qi>^T<^l5q) . (3.17) 

The isospin-one contribution can be represented in the form 

3 

-^NP,! = ^ 

^ i=l 

where 

Oi= {ql^T<^(y2.qWl^lTH) + {q'y^^-i5T<^azq){i-iyTH), (3.19) 

^>2 = {q'y^j5T'^o-3q){iji,j5TH) - {qjf^'y5T‘^a3q){t-fi,T°'t ), (3.20) 

^3= {q'y'^T'^(^3q)(J'yiM'y5T'^t) + {qi^i5T^cr3q){t-f^TH) . (3.21) 


2 ^' I,J=V,A 


cijOij 


(3.15) 
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Our isospin-zero and -one couplings c/j; Q are related to those of [22] by glcij/ml = cijj 
g'gCi/mf = < 7 ^ 02 / 771 ^ = d^Jb?^ glc^lmf = c'^^/A^. The relation to the coefficients 

of the operator basis (3.7) - (3.13) is given in [22]. To summarize, we use 

£nP = -CNP,g + -^^NP,q (3.22) 

The NP contributions to the coefficients of the tt spin density matrices (2.4) induced by 
interference with the tree-level QCD amplitudes of gg,qq —)• tt are listed in Appendix A. 

A side remark concerns bb —)• tt which receives, apart from the order QCD term, a SM 
contribution of order agCt from the interference of the t-channel kP-exchange - with possible 
anomalous Wtb couplings - and s-channel gluon exchange diagram. These SM terms, 
which are small, are taken into account in section 4. The contributions from anomalous 
Wtb couplings are neglected, because i) they are empirically known to be small (see the 
next subsection) and ii) the respective interference terms are suppressed by the bb parton 
luminosity. 

3.2 Top-quark decay 

The SM and NP contributions to the parton matrix elements that describe tt production 
and decay into dileptonic and lepton plus jets final states involve also the top-quark decay 
density matrices for t —)• Wb —)• ivb, qq'b. The top-decay vertex t —?■ Wb may also be 
affected by new physics interactions that can also be parametrized by anomalous couplings 
(cf., for instance, [20, 61]). We use in the following as spin-analyzers of the (anti)top quark 
only the charged lepton from decay, and we consider below only lepton angular 
correlations and distributions that are inclusive in the lepton energies. It is known [74-77] 
that these observables are not affected by anomalous couplings from top-quark decay if 
these couplings are small, i.e., if a linear approximation is justified. This is indeed the case 
in view of the present upper bounds on the moduli of these couplings that can be inferred 
from the measured IT-boson helicity fractions in top-quark decay [78, 79]. In other words, 
the observables which we analyze in the next section are affected only by possible new 
physics contributions to ti production, which we parametrize by the anomalous couplings 
discussed above. 

We use the results of [80] for the polarized semileptonic and non-leptonic top-quark decays 
in the SM at NLO QCD. 

3.3 Experimental bounds and domain of applicability 

The size of the anomalous couplings of the NP interactions described by (3.22) are con¬ 
strained by the experimental data on hadronic ti production which so far show no significant 
sign of new physics. A recent (re)analysis of the contributions of only the chromo-moments 
jit and dt to the ti cross section and confrontation with the Tevatron and LHC mea¬ 
surements yields correlated bounds. In terms of our (sign) conventions chosen in (3.4), 
Ref. [50] determined the allowed region (68% CL) in the jit, dt plane which may be rep¬ 
resented by the bounds —0.03 < jit < 0.02 and \dt\ < 0.15. Ref. [55] found a similar 
result: —0.046 < jit < 0.024 and \dt\ < 0.17 (95% CL). Ref. [68] took into account the 
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QCD corrections to the cross-section contribution of fit and obtained —0.019 < /if < 0.018 
(95% CL), putting dt = 0. The CMS experiment analyzed, for dileptonic tt events at 7 
TeV, the distribution of the difference of the leptonic azimuthal angles, and searched for 
a contribution of a non-zero fit- Using the predictions of [15] for interpreting their data, 
CMS obtained the bound [9] —0.043 < /if < 0.117 (95% CL). Tighter, but indirect bounds 
on dt and fit were obtained using the experimental upper bound on the neutron electric 
dipole moment [42] and data on rare B meson decays [42, 81], respectively. (Prospects for 
correlated bounds on the chromo moments in single top-quark production in the tW mode 
were recently discussed in [69].) 

A global fit to ti data using all dimension-six operators that are relevant to hadronic 
top quark production was recently performed in [70]. The operator set used in this reference 
contains also the operators contained in (3.22), albeit in a slightly different basis. One may 
use the correlated bounds of [70] to get order-of-magnitude bounds on the dimensionless 
couplings of the four-quark operators (3.14), to wit, |c 7 j|, |cj| < 0.4. 

In the next section we will compute the contributions of the NP interactions (3.22) to 
a number of observables, and we take into account in that section only the NP contribu¬ 
tions that are linear in the anomalous couplings. As the anomalous couplings correspond 
to non-renormalizable effective interactions, their contributions grow in the regime of large 
energy-momentum transfer to the ti system. In other words, the linear approximation can 
be applied with confidence only for relatively small anomalous couplings. In order to assess 
the domain of applicability of this approximation we have considered ti production at the 
LHC (14 TeV). We have computed, for the ti cross section and for a few spin correlation 
observables which will be introduced in the next section, the NP contributions to these 
observables which are linear and quadratic (respectively bilinear) in the anomalous cou¬ 
plings. Following Ref. [22] we define the range of applicability of the linear approximation 
by requiring that the contributions linear in the anomalous couplings are at least twice as 
large in magnitude as the bilinear ones. Application of this criterion yields \fLt\, \dt\ < 0.15 
and \cjj\ < 0.1. (We did not consider the contributions of the anomalous couplings Cj 
because they are subleading, see the next section.) Thus, as far as the chromo-moments 
are concerned, the empirical bounds discussed above are essentially within the domain of 
applicability of the linear approximation. 

4 Observables and results 

In this section we consider ti production at the LHC for center-of-mass energies 8, 13, 
and 14 TeV. First we determine the SM contributions and those of £np to the ti cross 
section ati and the ti charge asymmetry Ac- Then we introduce, for dileptonic and lepton 
plus jets final states, a set of correlation and polarization observables and compute their 
distributions. 

As far as the SM interactions are concerned we compute the four P-even and CP- 
even correlation observables that will be defined in section 4.2 at NLO QCD (order a^) 
[13] including the weak and mixed QCD-weak corrections of order a‘^,asa, 

[14, 15, 98-100]. We label these SM contributions by the acronym ‘NLOW’. The mixed 
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QCD-QED contributions to these spin correlation observables have not yet been calculated. 
We expect them to be small, of the order of a few percent of the QCD corrections. The 
combinations of polarization observables considered in section 4.2 are, apart from one com¬ 
bination, P-odd and/or CP-odd. These observables receive no QCD or QED contributions. 
Only one of these polarization observables, the sum of the transverse t and t polarization, 
receives, at the order of perturbation theory considered here, QCD contributions of order 
a^, and mixed QCD-weak and mixed QCD-QED contributions of order a^a, which are 
taken into account. (Contributions of order a^ag to the transverse polarization are very 
small and are not considered here.) The label ‘EW’ refers to the electroweak contributions. 
Our computational set-up is described in [14]. 

In the SM the total cross section is known to NNLO QCD (order af) including soft- 
gluon resummation at next-to-next-to-leading logarithmic (NNLL) order [72]. The mixed 
QCD-weak corrections to the cross section (NLOW) were computed in [98-102] and the 
mixed QCD-QED corrections were determined in [82]. The NLOW corrections to the LHC 
tt cross section are negative, about ~ —2% of the LO QCD cross section, while the QCD- 
QED corrections are positive, about ~ 1% according to [82]. The LHC charge asymmetry 
Ac defined below is computed at NLO QCD, including the mixed QCD-weak and mixed 
QCD-QED corrections [86]. These mixed EW corrections are known to be important for Ac- 
The computation of Ac to this order of perturbation theory, which is labeled ‘NLO+EW’, is 
present state-of-the-art. The Tevatron charge asymmetry was recently computed to NNLO 
accuracy [89]. 

Furthermore, we compute the interferences of the NP amplitudes induced by /Inp and 
the tree-level QCD amplitudes, but keep only the terms linear in the anomalous couplings. 
We employ the following input parameters: For the top-quark mass in the on-shell scheme 
we use mt = 173.34 GeV. Moreover, we use Pj = 1.3 GeV, mz = 91.2 GeV, mw = 
80.4 GeV, Pw = 2.09 GeV, mn = 125 GeV, a{mt) = 0.008, and we use the CTIO NLO 
parton distribution functions (PDF) [83]. This PDF set provides also the QCD coupling ag 
in the MS scheme. For the renormalization and factorization scale, we use fin = /Uj? = /x, 
fj, = m-t, and estimate scale uncertainties by varying /x between mt/2 and 2mt. 

4.1 Total tt cross section and charge asymmetry at the LHC 

Our results for are listed in table 2. The uncertainties refer to the scale uncertainties. In 
the linear approximation, only the P- and CP-even operators contained in £np contribute; 
the CMDM operator, whose contribution is largest, the isospin-zero operator Oyv^ and the 
P-even part of the isospin-one operator whose contribution is smallest. The smallness 
of the isospin-1 compared to the isospin-0 contribution to att, to the charge asymmetry 
defined below, and to the observables of the next section is due to the fact that a sizeable 
fraction of these contributions comes from the kinematic region near threshold. In this 
region, the contributions of uu —?■ tt and dd —)• tt are of the same order. The SM result 
for the cross section listed in table 2 is given here only for reference purposes. Comparing 
experimental results to predictions one should use the most precise available SM result on 
(Tff. This is the result of [72] mentioned above. Rather than performing a three-parameter 
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fit to cross section results at the Tevatron and the LHC (7 and 8 TeV) which is left to 
a future analysis, we make here a crude estimate of the bound on cyv by putting the 
anomalous couplings fit and ci to zero. Using, for 8 TeV, fiexp. = 241.5 ±8.5 pb [84] and the 
NNLO±NNLL QCD prediction [72, 85] fith = 253 ± 14 pb and neglecting the electroweak 
contributions, we get the 95 % C.L. bound —0.18 < cyv ± 0.09. The lower bound is outside 
of the condition \cvv\ ^ 0-1 for the linear approximation derived in section 3.3 if the other 
anomalous couplings are put to zero. 


Table 2: The SM and NP contributions to the tt cross section and the LHC charge 
asymmetry (4.1). 



NLOW 

oc cvv 

oc Cl 

oc fit 


8 TeV 

207.78;i|i 

246 

09 00 + 6.13 

OZ.oo_g 23 


cr[pb] 

13 TeV 

687.32;^8.| 

581.1611°^;! 

71 06 '’'^ 2-23 

‘ r-UD_i5 go 

2490.9ll5lf:?^ 


14 TeV 

8U.38+11-11 

654.65li^|20 

70 14+13.37 
<y-l+ 17.40 

2945.65is?|38 


NLO ± EW 

oc caa 

OC C2 


8 TeV 

(l.ll+[[;[!4) . 10-2 

U.4^D_o 017 

(9.26l°J2) . 10-2 

Tic 

13 TeV 

(0.75^^^) • 10-2 

u.o^i_o on 

(7.07l°:27) . 10-2 


14 TeV 

(0.66l°;°5) . 10-2 

n Q97+0.010 

U.OZ/_ooio 

(6.77l°:i) • 10“' 


Table 2 contains also our SM and NP prediction for the tt charge asymmetry at the 
LHC which is defined by 


^ cr{A\y\ > 0) - a{A\y\ < 0) 
^ (T(A|y| > 0) ± cj(A|y| < 0) 


(4.1) 


where A|y| = \yt\ — \yt\ is the difference of the moduli of the t and t rapidities in the 
laboratory frame. Our SM predictions, labeled ‘NLO ± EW’, contain apart from the mixed 
QCD-weak interaction corrections also the mixed QCD-QED corrections of order a^a [86]. 
We have computed Ac by expanding the ratio (4.1), see Eq. (4.22) below. The NLOW 
predictions in table 2 are in accord with the results of [86, 87] where different PDF sets were 
used. The NP contributions result from the isospin-zero operator Oaa and the P-even part 
of the isospin-one operator O^- These operators induce contributions to the differential 
cross section that are odd under interchange of the t and t momenta while those of the 
initial (anti)quark are kept fixed. Using the recent CMS measurement [88] with their 8 
TeV data, Ac = 0.001 ± 0.008, and putting the coupling C 2 to zero, we obtain the 95 % 
C.L. bound —0.06 < caa ± 0.01; cf. also [88, 93]. This result justifies a posteriori the use 
of the linear approximation. New physics contributions to the ti charge asymmetry ApB 
at the Tevatron were analyzed in terms of the above effective interactions in [22, 90-92]. 
The latest experimental results on the Tevatron asymmetry Aps are somewhat larger - 
but compatible - with the most precise SM results [89]. This implies caa > 0 and caa 
is bounded from above. The interaction described by Oaa is the effective field theory 
version of a number of NP models which were invoked to resolve the seeming discrepancy 
between the experimental results on the Tevatron App and the SM predictions. These 
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models include the axigluon model and flavor-changing t-channel Z' exchange in uu —)• tt] 
see [94] for a recent overview of the NP models discussed in the context of the tt charge 
asymmetry. 

The linear approximation for the contributions of the NP interactions to the cross sec¬ 
tion and the charge asymmetry is justified if \fit\ < |/it,max| — 0.15 and |c/j| < |c/j^max| — 
0.1, see section 3.3. We assume in the following that the latter bound holds also for 
anomalous couplings Cj. With these bounds we may make a crude estimate of the maximal 
difference between (future) experimental measurements of an observable and the SM pre¬ 
diction that will still allow the the use of the linear approximation. In the case of the ti 
cross section, we estimate this maximal difference by 

— \^tt Imax — loiCyy^jnaxI T |®2Cl,max| T lO's/tt^maxI j (4'^) 

where the coefficients a* are given in Table 2. (We take the values of the a* for ^ = mt.) The 
maximal difference 5Ac is determined in analogous fashion. In Table 3 our estimates of 
and 5Ac are given for the LHC (13 TeV). As these are crude estimates we have refrained 
from taking the scale uncertainties of the Oj into account. For the other LHC center-of- 
mass energies, the and 5Ac are obtained in analogous fashion using the coefficients of 

Table 2. 

Table 3: Moduli of maximal differences between future experimental result and SM pre¬ 
diction, for the LHC (13 TeV), for which the linear approximation can be applied. 


Observable 

no cut 

rriti < 500GeV 


440 

242 

6Ac 

4.1x10-2 

1.7x10-2 


Table 4 contains our predictions for and Ac for events with ti invariant mass < 
500 GeV. The mixed QCD-QED contributions to < 500 GeV), which are not 

included in this table, add ~ 1% to the cross-section values given in this table. The NNLO 
QCD corrections are not yet available. For Ac, the ratio of the isospin-one and isospin- 
zero NP contributions is approximately the same as in the inclusive case listed in table 2. 
In order to determine constraints on both anomalous couplings caa and C 2 from data on 
the LHC charge asymmetry alone, one may compute the differential charge asymmetry in 
narrower bins of or bins of the rapidity and transverse momentum of the ti system. 

In analogy to (4.2) one may estimate the maximal differences < 500 GeV) and 

8Ac{Mi^ < 500 GeV) that would allow the use of the linear approximation. They are given 
in Table 3. 

In the region where, for instance, the ti invariant mass is large, the linear approxi¬ 
mation in the NP couplings can no longer be justified just by assuming that the anomalous 
couplings are small, see section 3.3. While it is not difficult to compute ctii and Ac beyond 
this approximation, the following complication in the comparison with data arises: there is 
now a proliferation of contributions (being odd in the scattering angle y defined in section 2) 
which involve, apart from caa and C 2 , also the other anomalous couplings. Information on 
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both caa and C2 can be obtained either from considering Ac at different c.m. energies 
and the Tevatron asymmetry simultaneously, or from Ac and (future) measurements of the 
inclusive ti spin correlation and polarization observables defined in the next section. 


Table 4: The SM and NP contributions to the tt cross section and the LHC charge 
asymmetry (4.1) events with < 500 GeV. 


- 

NLOW 

oc cvv 

oc Cl 

oc /if 


8 TeV 

1 9Q VQ+lS-Sl 
izy. 


12.89l^:i 

480.35l}“.ui 

cr [pb] 

13 TeV 



23.981^61 

1455.30111^5 


14 TeV 


234 1 1+36.95 

26.121^88 

1702.521^^8.04 

- 

NLO EW 

OC caa 

oc C2 


8 TeV 

(o.geto;!]^) • 10-2 

u.zuy_Q Q44 

(4.38l[J:i) • 10-2 

Tic 

13 TeV 

(0.63;°:°^) • 10-2 

Q 14 .+0.007 
U--L+-L_o.007 

(2.89l°:{®) • 10-' 


14 TeV 

(o.5o;°:°^) • 10-2 

n 1 qo+0.006 

U.100_Q QQg 

{2.7ltlit) • 10-2 


4.2 Angular distributions and correlations for dilepton and jets events 

We analyze ti production at the LHC and subsequent decay into dileptonic final states, 

pp ^ t + i + X ^ + jets + Li™®, (4.3) 

and into lepton plus jets final states, 

pp ^t + i+ X ->f++ jets + , (4.4) 

pp ^ t + i+ X —jets + Li™'®®. (4.5) 


where f, I' = e, r. 

For the dileptonic events (4.3) we consider four-fold dilepton angular distributions. 
Without acceptance cuts, besides cuts on the ti invariant mass, these normalized distribu¬ 
tions have the form 


1 da 


1 

(47r)2 


(i + b; •!+ + B'2 • L 



(4.6) 


where dil = d cos 6d(f). The unit vectors £+, are the £~^ and directions of flight in 
the t and i rest frames, respectively, which are reached from the ti ZMF by rotation-free 
boosts. The coefficients, i.e., the vectors B'^, B 2 and the matrix C contain the information 
from the ti production and the t and i decay density matrices, which may be computed in 
the SM or in a NP model. As already mentioned in section 3.2, we consider the main SM 
top-decay mode t —)• LF6. We parametrize i+,i- in terms of spherical coordinates, but we 
allow for different polar axes, i.e., different reference axes a and b. We define 


cos 0+ 


cos 9- = i- ■ h . 


(4.7) 
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Integrating over the azimuthal angles we obtain the polar angle double distributions for a 
choice of reference axes a, b: 

—;- 7 ,—,- ^ = -fl + Bi cos 0+ + B 2 cos 0- — C cos 0+ cos 9- ). (4.8) 

a dcosO+dcosd- 4V / 


The sign in front of C is chosen such that in the absence of acceptance cuts [12, 13]: 


(^(a, b) = 


+ tT(U) - tj(ti) - 
cr(tt) + 


(4.9) 


where the ratio on the right-hand side of this formula is the double spin-asymmetry at the 
level of the t, i intermediate states. (The first (second) arrow refers to the spin state of the 
t (i) quark with respect to the axis a (b).) The top-spin analyzing power of the charged 
lepton from top decay is denoted by Ki. (In the convention used here, K£+ = ni- = K£.) At 
NLO QCD its value is K£ = 0.999 [80]. 

The polarization degrees of the ensembles of t and i quarks in ti events with respect to the 
reference axes a, b are given by 


P(a) = (2Si • a) , P(b) = (2Si • b), (4.10) 


where St, Sj denotes the t and t spin operator, respectively. They are related to Bi ^2 by 


.Bi(a) = P(a) K£ , B 2 (h) = -P{h) K£ . 


(4.11) 


That is, the relative signs in front of the B 12 in (4.8) are chosen such that in a CP-invariant 
theory and for the choice a = — b we have 


Bi = B2 . 


(4.12) 


Our choice of reference axes a and b which we use here is adapted to the orthonormal 
basis at the parton level introduced in section 2, but is not identical to it. As in section 2, 
we use the unit vector k which is the top quark direction of flight in the tt ZMF. Moreover, 
we use the direction of one of the proton beams in the laboratory frame, Pp, and define 
unit vectors fp and rip as follows: 


Pp = (0,0,1), fp = —(pp - ypk), 

rp 

?/p = Pp • k, 


% = —(PpXk), (4.13) 

rp 



Only in the case of 2 —)• 2 parton reactions and if the incoming parton 1 is parallel to Pp, 
the unit vectors defined in (2.5) are the same as those in (4.13). With the set (4.13) we 
define the reference axes a and b listed in table 5. The factors sign(yp) are required because 
of the Bose symmetry of the initial gg state. 

Instead of extracting the spin correlation coefficients C (for a choice of axes) from the 
double distributions (4.8) one may analyze the 1-dimensional distributions of the observ¬ 
ables 


^ = cos 0+ cos 0- . 


(4.14) 
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Labels are suppressed here, i.e., we consider in the following ^nn,^rr,^kk, ^rk,^rm etc. The 
one-dimensional distribution is obtained in standard fashion by inserting 


l = y d^5{^ — cos 9+cos 9-) (4-15) 

into the integral f d cos 9-^-dcos 9- over the double distribution (4.8). One arrives at 

which may be used for fits to unfolded experimental results. We recall that 

C = -9{0. (4.17) 

In this paper we apply the distributions (4.8) and (4.16) only to dileptonic final states. 
Integrating out one of the cosines in the double distribution (4.8) one gets 

= Vl-hfipa cos 6 i±y (4.18) 

a d cos Wj- 2 V / 

The slope of this distribution, Bi{B 2 ), measures the polarization of t (t) quarks with respect 
to the chosen reference axis. One may also compute/measure these polarizations using the 
relation Bi ^2 = 3 {cos9±). If no acceptance cuts are applied, the coefficient Bi {B 2 ) is the 
same both for the distribution of {i~) from dilepton and from lepton plus jets events. 

From now on we use the following notation. The label (a, b) refers to the choice of 
reference axes a and b from table 5. The correlation coefficient C{a, b) in (4.8) or (4.16) is 
associated with this choice of axes, and likewise i?i(a) and B 2 (b). 

In addition to (4.18), we compute also the top and antitop polarization with respect 
to another set of reference axes: 

B,{k*), B 2 {k*), Bi{r*), B 2 {r*), (4.19) 

where the labels k* and r* refer to the following vectors: 

k* : a = sign(A|y|) k , b = — sign(A|y|) k , (4.20) 

r* : a = sign(A|?/|) sign( 7 /p) fp , b = -sign(A|?/|) sign( 7 /p) fp . (4.21) 

Here A|y| = | 2 /t| — |yt| denotes the difference of the moduli of the t and t rapidities in 
the laboratory frame. The reason for introducing these additional polarization observables 
is that they are sensitive to the NP contributions from the operator Oav and from a P- 
odd combination of the operators Oj (defined in (3.16) and (3.19), respectively), while 
Bi^ 2 {k), Bi^ 2 {'k) project onto the contributions of Ova and O 3 . 

The coefficients C{a, b), Hi (a), and B 2 {b) (and sums and differences) are listed in 
table 6 . They can be measured with the 2-dimensional distributions (4.8) or with the 1- 
dimensional distributions of ^ab-, respectively with (4.18). Column 2 of this table indicates 
which coefficient of the ij —?• ttX production density matrices introduced in section 2 
contributes to the respective correlation coefficient. The type of interactions, i.e., the type 
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Table 5: Choice of reference axes. The unit vectors n^, fp and the variable yp are defined 
in (4.13). 


label 

a 

b 

transverse 

n 

sign( 2 /p) rip 

-sign(yp) rip 

r axis 

r 

sign(yp) fp 

-sign(yp) fp 

helicity 

k 

k 

-k 


Table 6: The correlation coefficients of (4.8), respectively of (4.17), (4.18) and (4.19), and 
sums and differences for different choices of reference axes. 


Correlation 


sensitive to 

C(n, n) 

^nn 

P-, CP-even 

C(r, r) 

^rr 

P-, CP-even 

C{k, k) 

^kk 

P-, CP-even 

C(r, k) + C{k, r) 

4 k 

P-, CP-even 

C{n, r) + C{r, n) 

^rn 

P-odd, CP-even, absorptive 

C{n, k) + C{k, n) 

4 n 

P-odd, CP-even, absorptive 

C(r, k) - C{k, r) 


P-even, CP-odd, absorptive 

C{n, r) — C{r, n) 

^k 

P-odd, CP-odd 

C{n, k) - C{k, n) 

-c^ 

P-odd, CP-odd 

Bi{n) + B2{n) 

1 ) 1 + _|_ i)^- 

P-, CP-even, absorptive 

Bi{n) - B2{n) 

— b^~ 

P-even, CP-odd 

Bi{r) + B2{r) 

bl^- + bl~ 

P-odd, CP-even 

B,{r)-B 2 {r) 

bl+ - bi- 

P-odd, CP-odd, absorptive 

B^{k) + B 2 {k) 


P-odd,CP-even 

Bi{k)-B 2 {k) 

4 ^ - 4- 

P-odd, CP-odd, absorptive 

Bi{k*) + B 2 {k*) 

4 ^ + 4 ~ 

P-odd,CP-even 

Bi{k*) - B 2 {k*) 

4 ^ - 4- 

P-odd, CP-odd, absorptive 

Bi{r*) + B2{r*) 

bl+ + bl- 

P-odd, CP-even 

Bi{r*) - B2{r*) 

bl+ - bl- 

P-odd, CP-odd, absorptive 


of contributions in the differential parton cross sections to which a correlation is sensitive, 
is given in column 3. 

Tables 7, 8 contain our results for the inclusive correlation coefficients C{a, b), Bi{a), B 2 {b) 
and Bi{a*), ^ 2 ( 6 *) for the LHC at 8 , 13, and 14 TeV, while the coefficients in tables 9, 10 
refer to events with tt invariant mass < 500 GeV. Several comments are in order. 

1. All normalized differential distributions, respectively the correlation coefficients are 
ratios. The P- and CP-even coefficients C are, in the SM at NLOW and to first order 
in the anomalous couplings, schematically of the form 

^ _ IVo + <5A^i + 5N^p 
gq + 5ai + (5(Tnp ’ 
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where Nq (<5A^i) and do ((5di) are the contributions at LO QCD (NLOW). A Taylor 
expansion at NLO in the SM and anomalous couplings yields 



(4.22) 


We use this expansion which is in the spirit of perturbation theory®. 

2 . As expected (cf. section 4.1), the tables 7-10 show that for all observables the 
isospin-1 NP contribution from Txp.q is about one order of magnitude smaller than 
the isospin-0 contribution. 

3. Within the Standard Model C{r,r) is, in the inclusive case, the smallest of the four 
P- and CP-even correlations. This can be qualitatively understood with simple tree- 
level arguments. At the LHC (8TeV) the t and t quarks in the inclusive tt sample 
are on average only moderately relativistic. At or very close to threshold the ti 
system produced at LO QCD by gluon fusion is in a state. Thus, not too far 
away from the threshold, the spin conhguration of the tt with respect to the r-axes 
hr and br (which point into opposite directions) is preferentially ff or J-i, while for 
relativistic tt the conhguration or dominates. This statement applies also to 
the gg —)> ti contribution to the helicity correlation C{k,k). The net effect for the 
inclusive sample is that gluon production contributes positively both to the correlation 
C(r, r) and C{k, k), but the correlation of the t and i spins with respect to the r-axes 
is signihcantly smaller than their correlation with respect to the helicity axes a = k, 
b = —k. The t, i produced by qq annihilation make a negative contribution both 
to C(r,r) and C{k, k), and it turns out that the magnitude of the contribution is in 
both cases approximately the same. This is because in this case, at LO QCD, ti near 
threshold is in a state. Thus both with respect to the r-axes and the helicity axes, 
the spin conhguration f), or \.'l dominates, both for non-relativistic and relativistic t, 
i. As a result, the contributions from gg and qq initial states almost cancel in the case 
of C{r,r), while in the case of C{k, k) the positive contribution from the gg channel 
is signihcantly larger in magnitude than the negative one from qq —)• ti. 

This line of arguments suggests that the ratio of C{r,r) and C{k, k) becomes larger 
if the ti sample is restricted to events not too far away from the ti threshold. This is 
corroborated by the numbers of table 9. The fact that the correlation C{n, n) is rather 
large in the SM can also be understood from total angular momentum conservation. 
The measurement of these four correlations would provide sufficient information to 
simultaneously determine or constrain the CMDM fit and the anomalous four-quark 
couplings cvv and ci. 

As an aside we remark that the NLOW and the CMDM [jit) contribution to the 
helicity correlation C(A;, k) at 8 TeV and the NLOW result for 14 TeV listed in table 7 
agree with those computed in [15] and in [14], respectively, with a different PDF set. 

® Ref. [68] analyzed the CMDM contribution /tt to the helicity correlation C{k,k) including the QCD 
corrections to this anomalous contribution. It was found that these QCD corrections do not improve the 
sensitivity to fit- 
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4. The opening angle distribution, a~^da/dcosip = (1 — D cosip)/2 [13], which is also a 
useful spin correlation observable, can be obtained from the above diagonal correlation 
coefficients. Because the reference axes defined in table 5 form orthonormal bases, we 
have [13] 

D = -l [C{n, n) + C{r, r) + C{k, k)] . (4.23) 

O 

This relation can be applied both to the SM and the NP contributions given in the 
tables below. The correlation coefficients D computed with (4.23) agree with those 
computed in [15] and [14] for 8 TeV and 14 TeV, respectively. 

5. One may also consider the spin correlations C{r*, k), C{k, r*), and C{k, k*), where the 
k* and r* axes are defined in (4.20) and (4.21). These observables are spin correlation 
analogues of the charge asymmetry Ac- They receive contributions from the entries 
of the ti production density matrices that are odd under interchange of the t and i. 
In the SM these correlations, which receive no contributions from the gg initial state 
and no contributions, at LO QCD, from qq initial states, are expected to be small. 
They receive NP contributions proportional to the anomalous couplings caa and C 2 - 

6 . As long as one considers CP-invariant interactions, C{r,k) = C{k,r). A non-zero 
difference C{r, k) — C{k,r) requires P-even but CP-odd absorptive contributions. In 
the SM at NLO there are no contributions of this type. Because our NP couplings 
are real and we do not take into account radiative corrections to the NP amplitudes, 
there are also no NP contributions to C{r, k) — C{k,r) in our approximation. 

7. The sums 

C{n,r) + C{r,n), C{n,k) + C{k,n) (4.24) 

probe P-odd, CP-even, naive T-odd absorptive (SM) contributions. The coefficients 
C{n,r) and C{n,k) receive SM contributions from P-odd absorptive parts of the 
mixed QCD-weak corrections. We have computed the respective 1-Ioop absorptive 
contributions to the gg —)• tt and qq —)• tt spin density matrices. The resulting effects 
on (4.24) are very small, cf. tables 7 and 9. Because all anomalous couplings are 
real in the effective Lagrangian approach and because we do not take into account 
radiative corrections to the NP amplitudes, there are no NP contributions to (4.24). 

8 . The SM contributions referred to in the previous item cancel in 

C{n,r) — C{r,n), C{n,k) — C{k,n). (4.25) 

These differences probe P-and CP-odd (dispersive) NP contributions, that is, non¬ 
standard CP violation. These differences in the correlation coefficients can be related 
to CP-odd triple correlations. Defining 

= (£+ xL)-k, or = sign(yp)(4 x L) • fp, (4.26) 

the following relations hold: 

C{n,r)-Cir,n) = 9{Of^), C{n,k) - C{k,n) =. (4.27) 
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Alternatively one may consider associated CP asymmetries, which are easier to mea¬ 
sure. They are related to the expectation values of (4.26) by [15, 95]: 


CP _ Nu{OfP>Q)-Nu{0?^ <^) _ , 


(4.28) 


Formulae (4.27) and (4.28) hold if no cuts on the dileptonic final states other than a 
cut on the ti invariant mass are applied. The respective numbers in Tables 7 and 9 
show that C{n, r) — C{r, n) is the most sensitive spin correlation observable to search 
for a top-quark chromo-electric dipole moment dt- The CP asymmetry is less 
sensitive to the CP-violating anomalous couplings than the left-hand sides of (4.27), 
but experimentally more robust. The measurement of the two correlations (4.27) or 
the two asymmetries (4.28) would allow to simultaneously determine or constrain the 

CEDM dt and the anomalous CP-odd gluon-top-quark coupling C(_). 

Putting the coupling C(__) = 0, the expectation value at 8 TeV computed with 

the 8 TeV value for C{n, r) — C{r, n) given in table 7 agrees with the result of [15]. 


9. As already emphasized above, for pp collisions a classification of observables with 
respect to CP is, strictly speaking, not possible because \pp) is not a CP eigenstate. 
However, there will be no SM contributions from qq and gg initial states to the CP- 
odd observables listed in the previous item and in (4.30) - assuming no acceptance 
cuts or cuts which are CP-invariant. Moreover, within our approximation we found 
no contribution from P-violating but CP-even terms induced by gq production. 


10. As to the longitudinal top-polarization coefficients B{k), B(r): The sums 


Bi{k)+B 2 {k), Hi(r) + H 2 (r), Bi{k*) + B 2 {k*), B^{r*) + B 2 {r*) (4.29) 


receive contributions from the P-odd weak-interaction corrections, that is, from the P- 
odd terms of the NLOW corrections. However, the SM-induced t and i polarizations 
are very small to tiny, as the numbers in table 8 and 10 show. (The SM result 
for Bi[k) + B 2 {k) at 8 TeV listed in table 8 agrees with the result computed in 
[15].) The sums (4.29) receive also contributions from the P-odd NP four-quark 
operators. Thus, the measurement of the four polarization observables (4.29) would 
allow to simultaneously determine or constrain four anomalous four-quark couplings, 
for instance cyA-, cav-, C 2 , and £ 3 , assuming that input on ci is provided from the 
measurement of the four CP-even spin correlations, see item 3 above. 

In the differences 


B,{k)-B 2 {k), Hi(r)-H 2 (r), Bi{k*) - B 2 {k*), Bi{r*) - B 2 {r*) , (4.30) 

which are P- and CP-odd and Tjv-even, the SM contributions to Bi ^2 cancel. Non¬ 
zero differences require, in addition to CP-violating NP contributions, also absorptive 
parts. Because our anomalous couplings are real and we do not take into account 
radiative corrections to the NP contributions, these differences remain zero within our 
approach. In [15] the chromo moments were interpreted as form factors which may 
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have an absorptive part, and the differences of the longitudinal t and t polarizations 
were computed in terms of Imdt- 

11 . As to the combinations 

Bi{n) ±B 2 {n). 

The sum is generated by P- and CP-even absorptive terms in the differential cross 
section. In the SM such terms arise from a number of one-loop amplitudes for gg, qq —)• 
ti which involve parton, weak gauge-boson, Higgs boson, and photon exchange. The 
QCD contributions given in tables 8 are in accord with the 8 TeV result of [15]. (In 
this work, the vector perpendicular to the scattering plane was used in un-normalized 
form and a different PDF set was employed.) In addition we list in tables 8 and 10 also 
our results for the electroweak contributions, which arise from the absorptive parts 
of the mixed QCD-weak interaction and QCD-QED contributions to gg, qq ti. 
They turn out to be non-negligible in comparison to the QCD contribution. Thus, for 
the inclusive ti sample, the transverse polarization of the ensemble of t and i quarks 
with respect to the reference axis sign(yp)np is Pnt = Pnt — 0.5%. In the difference 
Bi(n) — B 2 {n) the SM contributions cancel. The difference is sensitive to P-even, 
but CP-odd dispersive NP contributions, which in the framework of our effective 

Lagrangian approach are generated at Born level by the operator _defined in 

(3.5). 

12. Finally we add a remark on the scale dependence of our SM predictions for the 

correlation and polarization observables. As the numbers given in Tables 8 and 10 
show, the parity-odd combinations of the B coefficients are very small in the SM but 
exhibit large scale variations relative to the values for g = mt- The reason for this is 
as follows. For the numerators N of the B coefficients that are generated by the P-odd 
weak interaction contributions, one gets N(g = 2mt) > N[g = mt) > = mtjT), 

due to factorization-scale dependent logarithms in some of the qg induced mixed QCD- 
weak corrections [100], while for the denominators D of the B coefficients, which are 
taken into account at lowest order QCD, one has D(g = 2mt) < D{g = mt) < 
D{g = mtj2). This produces the large scale variations in the ratios B = N/D. 
The P- and CP-even correlation observables given in Tables 7 and 9 are generated 
predominantly by QCD. Here the scale variations are moderate to small and, as is 
well known, the Taylor expansion (4.22) reduces the scale dependence of the results at 
NLO QCD with respect to LO. The weak interaction contributions with their rather 
large scale dependence do not spoil this because they are small compared to the 
QCD contributions. The correlation C{r,r) is special in that there are cancellations 
among separate contributions, keeping C{r, r) rather small and more sensitive to scale 
variations than C{n,n) and C{k,k). 

The above spin-correlation and polarization observables probe all entries of the top- 
spin dependent parts of the ti production density matrices. In addition, the various NP 
contributions are disentangled to a large extent, as the above discussion and tables show. 
Moreover, a measurement of this set of spin observables - and experimental results on 
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Ac, which yields information on the coupling caa ^ would provide sufficient input for 
determining all anomalous couplings of the effective Lagrangian (3.22). As mentioned 
above, the spin correlations C{r*,k), C{k,r*), and C{k,k*) are also sensitive to caa- 

An important issue is the range of validity of the linear approximation for the NP 
contributions to the correlation and polarization observables. In analogy to the procedure 
we applied to the tt cross section and the charge asymmetry Ac in section 4.1, - cf. Eq. 
(4.2) and Table 3 - we estimate the maximal difference between (future) experimental 
measurements and SM prediction of an observable by adding the moduli of the products of 
the respective contributions given in Tables 7-10 and the maximal values of the anomalous 
couplings, derived in section 3.3, for which the linear approximation can be applied. The 
results of these estimates are given, for the LHC (13 TeV) in Table 11. The estimates for 
8 and 14 TeV can be obtained in completely analogous fashion. The numbers of Table 11 
show that the range of validity of the linear approximation is quite large. 

One may compute the above spin observables also for tt events with large tt invariant 
mass, for instance M^i > 500 GeV. While there is no computational difficulty, one should, 
as discussed above, in this case abandon the linear approximation in the anomalous cou¬ 
plings. This leads to a proliferation of NP contributions to each of the spin-correlation and 
polarization coefficients, which will make a comparison with (future) data rather intrans¬ 
parent. One may pursue a pragmatic approach, namely, analyze only selected observables 
beyond the linear approximation, in particular those which turn out to be measurable with 
moderate or small experimental uncertainties. Such an analysis is beyond the scope of this 
paper and will be left to a future investigation. 
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Table 7 : The correlation coefficients C at NLOW in the SM and the non-zero contributions 
of the dimensionless anomalous couplings defined in section 3. The central values correspond 
to the scales fi = mt- The uncertainties refer to fj, = mt/2 and /r = 2mt. 


- 

NLOW 

oc cv-y 

oc Cl 

oc fit 

C{n, n) 

8 TeV 

13 TeV 

14 TeV 

n QQ9+0-002 
U.OOZ_Q QQ2 

U.OZU_Q gQ2 
n Q9 c;+0-002 

-0.150l[j;;]J“ 
(-7.99l°|°) • 10-2 

(-7.21+°f) . ^0-2 

(-1.76l°13) . 10-2 
(-8.211°;^^) • 10-3 
(-7.201°!^) • 10-3 

1 8*^1 +0.038 

i.001_Q Q3g 

9 n9c;+0.025 

z.uzo_Q Q24 

2 047''“°'°22 
Z.U4r-o 022 

C{r, r) 

8 TeV 

13 TeV 

14 TeV 

U.UOO_0 007 

-'■-0.006 
U.U/Z-o 006 

— 0 Q‘14+'^ ''33 

u.ao^_0 037 

—0 697'’'^'*^27 
u.uy (_o 027 

U.DDO o 020 

n 1 99+0.006 

U.1ZZ_Q QQg 

(-8.461°;") • 10-2 
(-7.99l°:39) . 10-2 

9 q-| q+0.032 
^•oiO_0 034 

9 47^+0.020 
Z.4/0-0 019 

^•ciyo 0 049 

C{k, k) 

8 TeV 

13 TeV 

14 TeV 

n *^1 O + 0.003 

U.OiO_Q QQ2 

n ‘^‘^1 + 0-002 

U.001_Q QQ2 

n ‘^‘^1 + 0.002 

U.001_Q QQ2 

i.om_0 053 

1 o-i q+0.039 
r-^ro_o 040 

— 1 169'’'^3537 
r-roy o 03g 

-0.2091°:°°° 

U. 10 i_o 007 

-0.1431°:°°° 

U.yD 4_0 009 

0 qi7+0.006 
'^•'^■'-'-0.006 

u.yii-o 005 

C(r, k) + C{k, r) 

8 TeV 

13 TeV 

14 TeV 

—0 226^°'°°^ 
U.ZZU _0 004 

U.ZUD _0 002 

U.ZU ^-0 002 

—0 430'’'^-^36 
U.^OU_o 04g 

— 0 306'''^3534 
U.OUD_o 044 

—0 290'’'^'°32 

u.zyu o 042 

(-5.42l°;|) • 10-2 
(-3.58l°;2°) . 10-2 
(-3.34l°;33) • 10-2 

‘ ^'J-O.OIO 

n 740+0.001 

O'. ‘"'^U- 0.002 

0 74 Q+ 0.002 

U. /4U-0 004 

- 

NLOW 

C{n, r) + C{r, n) 

8 TeV 

13 TeV 

14 TeV 

(1.03l°;°l) • 10 -=* 

(1.06l°;OJ) • 10-3 
(1.07l°;°3). 10-3 

C(n, fc) + C(/c, n) 

8 TeV 

13 TeV 

14 TeV 

(2.32l°:°J) • 10-3 
(2.15+0°^) • 10-3 
(2.12l°;03) . 10-3 

- 

(X C() 

oc dt 

C(n, r) — C(r, n) 

8 TeV 

13 TeV 

14 TeV 

1 *^91 + 0.021 

1.0Zi_Q Q2Q 

1 99^^+0.007 

l.ZZU_Q QQ 4 

1 91 '7+0.005 
■‘-•^-‘-^- 0.002 

q 7r\i:-l-0.084 

O. (UO -0 087 

—4 14Q+0.053 

4:.140-0 056 

—4 ;intc+0.052 

"i.ryo 0.052 

C(n, k) — C{k^ n) 

8 TeV 

13 TeV 

14 TeV 

J--ODO_0.058 

9 1 ^7+0.108 
^•-*-^^-0.087 

9 9 Qf^+ 0.112 

Z.ZOU_Q Q9;L 

n 791 + 0.021 
' ^-'--0.019 

-0.8001°;°°° 

-0.8051°;°°^ 
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Table 8: The correlation coefficients B at NLOW in the SM and the non-zero contributions 
of the dimensionless anomalous couplings defined in section 3. The central values correspond 
to the scales fi = mt- The uncertainties refer to fj, = mt/2 and /r = 2mt. 


- 

NLOW 

OC cvA 

OC C3 

Bi{r) + B 2 {r) 

8 TeV 

13 TeV 

14 TeV 

(2.9l};^) • 10-3 
(3.212-3) . 10-3 
(3.312-3) . 10-3 

U.ziu_o 009 

0 iqq+O-008 

u.iyy o,oo8 

(3.89l°;24) . 10-2 
(2.55l°:}l) • 10-^ 
(2.39l°;}3) . 10-2 

Bi{k) + B2{k) 

8 TeV 

13 TeV 

14 TeV 

(6.111^) • 10-3 
(8.0ll^) • 10-3 

(8.310) • 10"' 

^•rUU_o 078 

1 607'’'3-031 
i.DUr_o 052 

1 p;42+0-047 

l.O^Z 0 047 

0.283lO;°|t 

0.2011°;°“ 

0 iqi+0-008 

u.iai o 008 

- 

NLOW 

OC CAV 

OC Cl — C2 -|- C3 

Bi{r*) + B 2 {r*) 

8 TeV 

13 TeV 

14 TeV 

< 10-3 

< 10-3 

< 10-3 

r.iDD_o 060 

^■'^^-0.039 

0 746+3-037 

o’-'+O-o.ose 

0 244'’'°-°3^ 

U.Z^4:_o 014 

0 164+°'°°° 

U.iD4:_o 009 
n 1 I14+0-008 

U.i04:o 008 

Bi{k*) + B2{k*) 

8 TeV 

13 TeV 

14 TeV 

< 10-3 

< 10-3 

< 10-3 

1 204'’'3034 

i.ZU4_0 052 

n qoo+0.034 
U.05<3_q Q 34 

n Q/19+0-032 

U.042 q Q32 

0 211 +0-013 
U.Z01_o 012 

n 1 qo+ 0.008 
U.lOO_o 008 

n 174+0-007 
U-l'^-O-OOT 

- 

NLO QCD 

EW 

Bi{n) + B 2 {n) 

8 TeV 

13 TeV 

14 TeV 

( 6 . 93 IO 2 ) • 10-' 
(7.45lO^) • 10-3 
(7.52lO^) • 10-3 

( 3 . 44 III;?) . 10-3 

(3.891^6) . 10-3 

(3.941°;“) . 10-3 

- 

OC C(+) 

Bi{n) - B 2 {n) 

8 TeV 

13 TeV 

14 TeV 

4.148003 

4.87600^ 

4.97500^ 
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Table 9: The correlation coefficients C at NLOW in the SM and the non-zero contributions 
of the dimensionless anomalous couplings defined in section 3 for < 500 GeV. 


- 

NLOW 

OC cvv 

OC Cl 

OC p,t 


8 TeV 

u.oyu_Q gQ2 

n *^1 7+0.012 
^•'^-‘-'-0.013 

(-3.89;0;20) . 10-2 

i-ob4_Q Q34 

C{n, n) 

13 TeV 

0 4nQ+0-002 

U.'^UO_Q 002 

U.ZOU_o 040 

(-2.59;°!^) • 10-2 

2 n4i+0'*^24 

Z.U+i-o 023 


14 TeV 

0 4QQ+0-002 

U.4UO 0 002 

u.ziy o oog 

(-2.43;o13) . 10-2 

2 062'’'®3)^3 
^•UUZ o 022 


8 TeV 

n 1 Q7+0.007 
'-^•-‘-’^'-0.005 

n S91 +0.030 
U.OZi_Q Q32 

-0.10+liZ 

^•OUU_0 028 

C{r, r) 

13 TeV 

U.lO5_0_003 

—0 5qQ+°-02‘^ 

u.ayy_o 025 

(-6.78;0;i) • 10-2 

^•oru_0 020 


14 TeV 

n 1 71 

'-^•-'-'-'■-0.003 

n ^71 +0.022 
-'--0.023 

(-6.36;0;33) • 10-2 

9 c;97+0'019 
^•'-’^'-0.018 


8 TeV 

n 41Q+0-003 

U.'^iO-O 003 

r-UDU_o 040 

f) 1 01 +0.006 

U.lOl_o 007 

o.928i[j:;]3" 

C{k, k) 

13 TeV 

n AAf^+O-007 
U.44D_o 008 

n 777+0.030 
' ^-0.031 

(-8.80;°!^) • 10-2 

0.863l[J:°?9 


14 TeV 

U.448 0 008 

—0 741+0-028 
'^■'--0.029 

(-8.28;o13) • 10-2 

n cc:c;+0.008 

U.OOO 0.010 


8 TeV 

U.ZOU_o 004 

—0 430+°-*^^® 

U.^OU_o 040 

(-5.30;0;i) • 10-2 

0 980'''° '9i4 
U.y5U_o 014 

C(r, k) + C{k, r) 

13 TeV 

n 91 9+0-010 

U.Z1Z_Q QQ2 

—0 41R+0-012 
U.OiD_o 042 

(-3.56;°i3) • 10-2 

1 n54'’'^'°3o 
i-UO4_0 010 


14 TeV 

—0 210^®'^®^ 
U.ZiU_o 002 

U.OUU o 042 

(-3.36;oi6) . 10-2 

r-UDZ o 010 

- 

N 

LOW 


8 TeV 


< 10-3 


C(n, r) + C{r, n) 

13 TeV 


< 10-3 



14 TeV 


< 10-3 



8 TeV 


(2.64l0;[]^) • 10-3 


C(n, A;) + C(/c, n) 

13 TeV 


(2.87l0;[]^) • 10-3 



14 TeV 


(2.9110;°2) . 10-3 


- 

OC C() 

OC dt 


8 TeV 

—0 f)R7+0-035 

U.DD^-o 035 

—3 74f)+0'*^3® 

'J-‘^'3_oog4 

C(n, r) — C{r, n) 

13 TeV 

-0.481" 

4.ZU5_q Qg2 


14 TeV 

-0.459l°;[]i 

-4 261+°'°37 

+.ZDi o 059 


8 TeV 

0.719 

+0.036 

-0.035 

— 1 oq8+®*532 

r-uyo_o 034 

C(n, k) — C{k, n) 

13 TeV 

0.908 

+0.022 

-0.022 

1 97^^+0.023 
l.z/^-0 023 


14 TeV 

0.929 

+0.020 

-0.020 

i.zyct-o 022 
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Table 10: The correlation coefficients B at NLOW in the SM and the non-zero con¬ 
tributions of the dimensionless anomalous couplings defined in section 3 for < 500 
GeV. 


- 

NLOW 

OC CVA 

OC C 3 

Bi{r) + B 2 {r) 

8 TeV 

13 TeV 

14 TeV 

(2.5l|?) • 10-3 
(2.9l?;3) • 10-3 
(2.9+1?) • 10-3 

0 196+*^'°30 
u.iyD_o 010 

0 1 Q 4+0006 

U.i0'4_o 006 

n 1 o7“I"0-006 
'^•-'-^^-0.006 

(2.43l0;|0) • 10-2 
(1.52l0;00) . 10-2 
(1.42l0;03) • 10-2 

Bi{k) + B2{k) 

8 TeV 

13 TeV 

14 TeV 

( 2 . 6 l?;^) • 10-3 
(3.5l|^) • 10-3 
(3.6l|?) • 10-3 

n 070+0 050 
U.y (y_0 049 

0 669+3-^32 

U.DDy_o 031 

0 632+3-°^® 

U-UOZ o 029 

0 191 +0.008 
u-r^r_ 0.008 

(7.6010;^3) . 10-2 

(7.09lO;^?) • 10-2 

- 

NLOW 

OC CAv 

OC Cl — C 2 -|- C 3 

Bi{r*) + B 2 {r*) 

8 TeV 

13 TeV 

14 TeV 

< 10-3 

< 10-3 

< 10-3 

0.988l+[;i 

0.65910:034 

0 620+0-^31 
U-OZU o 031 

0 207+0-042 
4’-^'91_o 012 

fl 1 or+0.008 

U-roo_Q QQ 7 

0 1 26+0-^407 

U.iZO Q QQ 7 

Bi{k*) + B2{k*) 

8 TeV 

13 TeV 

14 TeV 

< 10-3 

< 10-3 

< 10-3 

0 7fitC+0-U42 
U. /DO_o 041 

0 51 1 +0-026 
U.Oii_o 025 

n zLRI +0-024 

U.45l o o24 

0 160+0-099 

u-rDU_Q Q^g 

0-10510:000 
(9.8310:00) . 10-2 

- 

NLO QCD 

EW 

Bi{n) + B 2 {n) 

8 TeV 

13 TeV 

14 TeV 

(4.121+t^) . 10-3 
(4.16l°+5^ • 10-3 
(4.16l°+5^ • 10-3 

(3.1410-09) . 10-3 

(3.501^3^) • 10-3 
(3.551^1) • 10-3 

- 

OC C(+) 

Bi{n) - B 2 {n) 

8 TeV 

13 TeV 

14 TeV 

^■"±y»_0.051 

9 7c:/(+0.036 
‘^^-0.038 

9 7cp:+0-034 

z. /Od_o 036 
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Table 11: Moduli of maximal differences between future experimental results and SM 
predictions, for the correlation and polarization observables at the LHC (13 TeV), for which 
the linear approximation can be applied. 


Observable 

no cut 

mil < 500GeV 

5C{n, n) 

0.31 

0.33 

6C{r, r) 

0.44 

0.44 

5C{k, k) 

0.27 

0.22 

S[Cir, k) + C{k, r)] 

0.15 

0.19 

S[C{n, r) — C{r, n)] 

0.74 

0.68 

5[C{n, k) - C{k, n)] 

0.34 

0.28 

5[B^{r) + B 2 {r)] 

2.4x10“2 

1.5x10“2 

5[Bi{k) + B2{k)] 

0.18 

7.5x10-2 

5[Bi{r*) + B 2 {r*)] 

0.13 

0.11 

5[B^{k*) + B2{k*)] 

0.14 

8.3x10-2 

5[Bi{n) - B 2 {n)] 

0.49 

0.28 


5 Summary and conclusions 

We have defined a set of correlation and polarization observables with which hadronic tt 
production dynamics can be probed in detail, that is, all entries of the top-spin depen¬ 
dent parts of the hadronic ti production density matrices, in dileptonic and semileptonic tt 
events at the LHC. Our observables allow, in particular, to disentangle contributions from 
interactions which conserve/violate parity and/or CP. We have computed these observables 
within the Standard Model at NLO QCD including the mixed QCD-(electro)weak contri¬ 
butions. In addition we have analyzed possible new physics contributions in terms of an 
effective Lagrangian that contains all gauge-invariant effective interactions of mass dimen¬ 
sion < 6 being relevant for the parton processes at hand. Existing constraints on anomalous 
couplings which describe the strength of these effective NP interactions justify that only 
NP terms linear in the anomalous couplings are taken into account in the computation of 
inclusive (spin) observables for the LHC. Irrespective of these constraints we have estimated 
the range of validity of the linear approximation for the observables considered in this pa¬ 
per. Some of our polarization and spin correlation observables allow for direct searches of 
non-standard P- and CP violation in tt events and will provide, once measurements will 
have been made, direct information on respective interactions which may affect hadronic tt 
production. The measurement of the complete set of our inclusive spin observables, includ¬ 
ing Ac, would allow to determine, respectively constrain all anomalous couplings contained 
in Tnp- Instead of Ac, one may consider the observables C{r*,k), C{k,r*), and C{k,k*) 
which are the spin correlation analogues of the charge asymmetry. 

Eventually, it will also be of interest to compute and measure the above tt spin cor¬ 
relations and t and t polarizations in the kinematic regime where the energy-momentum 
transfer to the tt system is large, for instance for tt events with > 500 GeV. In this 
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regime the contributions of the effective NP interactions grow in general, due to their non- 
renormalizable nature. In this kinematic regime one should abandon the linear approxima¬ 
tion, but then each observable receives a large number of contributions from the anomalous 
interactions contained in £np- Such an analysis should be done at some point for selected 
observables, once experimental feedback on measurement precisions will be available. 


A NP contributions to the tt spin density matrices 


The tt production spin density matrices to the parton reactions gg, qq —)• tt{g) and gq{q) —)• 
tiq{q) were computed at NLO QCD (dispersive contributions) in [12, 13], the contributions 
of the NLO QCD absorptive parts were determined in [96, 97] and the contributions of the 
mixed QCD-weak corrections (dispersive contributions) in [98-100]. The spin-independent 
term of the mixed QCD-weak corrections to qq,gg —?■ ti, i.e., the respective contributions to 
the coefficient functions of eq. (2.4), were computed also in [101, 102] and, without the 
infrared-divergent box contributions to qq —?■ ti and the corresponding real gluon radiation, 
in [103], and these results agree with ours. In addition, we have computed the contributions 
of the absorptive parts of the mixed QCD electroweak corrections. Here we do not list 
the latter contributions, for the sake of brevity and because they induce only very small 
effects (see section 4.2). Semileptonic and non-leptonic decays of polarized top-quarks are 
incorporated at NLO QCD [80]. 

Here we list the contributions to the tt spin density matrices (2.4) in the decomposi¬ 
tion (2.6), (2.7) of the NP interactions described by the effective Lagrangian (3.22). As 
emphasized above, we take into account only the interferences with the tree-level QCD 
amplitudes of gg, qq —)• tt, i.e., only terms linear in the anomalous couplings. We recall that 
all anomalous couplings are real dimensionless parameters. 

In the following, s = {pi + P 2 )‘^, where pi ^2 are the 4-momenta of the initial partons. 
Furthermore we define (cf. (2.5)) 


2; = 


2mt 


/3 = Vl - 2;2 , 


y = pk. 


(AT) 


A.l NP contributions to gg —)■ tt 

For the non-vanishing NP contributions to the coefficients of the gg —?• tt spin density 
matrix defined in (2.6), (2.7) we get in the linear approximation, dropping the superscript 

I = gg- 
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A.2 NP contributions to qq tt 

For the non-vanishing NP contributions to the coefficients of the qq —>■ tt spin density matrix 
defined in (2.6), (2.7) we get in the linear approximation, dropping the superscript I = qq: 
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The terms proportional to At) dt, and C(_) apply to all q ^ t. The terms proportional to 

the couplings cx apply only to the first quark generation; the upper (lower) sign refers to 
uu {dd) in the initial state. 
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